Effects of non-denumerable fixed points in finite dynamical systems 
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The motion of a spinning football brings forth the possible existence of a whole class of finite 
dynamical systems where there may be non-denumerably infinite number of fixed points. They defy 
the very traditional meaning of the fixed point that a point on the fixed point in the phase space 
should remain there forever, for, a fixed point can evolve as well! Under such considerations one 
can argue that a free-kicked football should be non-chaotic. 

PACS numbers: 05.45.-a,47.52.+j,01.80.+b 



I. INTRODUCTION 

Any elementary standard textbook [l| on non- linear dy- 
namics would show that for a finite dynamical system, 
may be non-linear, one always gets a set of finite number 
of fixed points. Interestingly, a seemingly overlooked case 
in the existing literature had been the one in which the 
parameters in such a system could depend on the initial 
conditions; some very recent highly rigourous mathemat- 
ical works (see Q and the references therein) on systems 
that contain a hidden parameter and thus, continuous 
one-parameter families of equilibria, have been done. As 
we shall show in this paper, such a finite dynamical sys- 
tem is of real interest, for, it shows an array of weird 
properties such as evolution of the fixed points, existence 
set of non-denumerable fixed points etc. Moreover, as we 
shall see, such systems are not just of non-realistic ped- 
agogical and academic interest but these systems can in 
fact exist in reality e.g., in the case of a sphere projected 
with a spin i.e., say a football or a tennis ball hit to im- 
part both the linear and angular velocities on them. 
As we know that the archy path of a kicked football is 
non-chaotic. The model equation used to describe it, 
as we shall see, is in fact such a dynamical system. As 
an application of these types of dynamical systems, it 
will be shown using numerics that the very existence of a 
bounded set of non-denumerably infinite fixed points aris- 
ing in the case of the model equation causes the path of 
the football to be predictable or non-chaotic. Thus, the 
natural existence of peculiar dynamical systems with pa- 
rameters depending on initial conditions should encour- 
age researchers to explore such systems both theoretically 
and numerically. 

In the section (II) , we shall recall the true reason for the 
bending of a ball during its flight and show how the wrong 
reason (Bernoulli's principle) gives a quantitative estima- 
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tion of the bending using which we frame the model equa- 
tion in section (III) that may even model the effect of the 
reverse Magnus effect. The model equation is obviously 
a set of three 1st order non-linear coupled differential 
equations with a parameter depending on initial condi- 
tions, an issue which is the main point being studied in 
this paper. In section (IV), we study the model equation 
analytically and numerically to show that the football's 
motion is non-chaotic owing to certain peculiarities of the 
model equation. The issue of striking peculiarities in the 
model equation is further taken up in the section (V). 



II. WHY FOOTBALL SWERVES? 

We all have seen expert footballers making weird 
swerving free-kicks by imparting spin on the football, 
little knowing, even qualitatively, how exactly it occurs. 
Though the bending of the path of a rotating sphere mov- 
ing through a fluid {e.g. air) is often attributed to the 
so-called Magnus-effect Q (or Robins-Magnus effect [J], to 
be historically unerring) owing to the Bernoulli's princi- 
ple, the process causing the bending of the spinning ball 
is very much complex even qualitatively, let alone the 
dream of an exact quantitative equation. 
How does Bernoulli's principle [lj, [3 explain the bend- 
ing of a rotating ball? To clarify how exactly the logic 
goes is given below with the help of an anti-clock-wisely 
rotating sphere of radius a moving forward in the direc- 
tion with uniform velocity v transverse to the direction 
of the axis of rotation, which passes through the center 
of the sphere, in a fluid (say, air of density p) at rest. 
(We call it 1-D case, for obvious reason.) The angular 
velocity w.r.t to the center be ui. Of course if one fixes 
oneself at the center of the sphere, then in that frame the 
sphere is rotating and the air is sweeping past the ball 
uniformly with velocity v (Fig. 1). The fluid at point P 
on the streamline adjacent to the periphery of a sliced 
disc (of radius r and infinitesimal thickness da) of the 
sphere is supposed to have velocity v of the air in the x- 
direction and — LursmO and — Lur cos in the x-direction 



2 




Front-view 



FIG. 1: Flow past a rotating sphere. Here, side-view of 
an infinitesimal thick disk or radius r has been shown with 
the adjacent streamlines (shown exaggeratedly away from the 
body). The front-view shows the necessary geometry needed 
to carry out the integration. 

and y-direction respectively due to the rotation of the 
sphere. The application of Bernoulli's principle suggests 
that pressure p at point P is: 

p = pcuvr sin 8 — — pui 2 r 2 (1) 

The total force F (Magnus force) on the sphere in the -y- 
direction may be obtained by integrating p(rd9 sin 6) (da) 
over the upper and the lower streamlines (remembering 

r= vV-O): 

/•cr— +a 1*9— 2n j 

F = I / prsm8d8da = -ira 3 pujv (2) 

J <r=-a Je=o 3 

Although seemingly okay, any careful reader must have 
noted that there are flaws in the argument! Bernoulli's 
principle, in the form used here, is applicable for the 
steady, incompressible and inviscid fluid; the first two re- 
striction can be justified but the last one can't be justified 
near the boundary. One may note, the no-slip condition 
has been implicitly assumed in the derivation by allow- 
ing the fluid adjacent to the boundary to move with the 
sphere but this is in direct opposition to the concept of 
ideal fluid in which Bernoulli's principle is applicable! 
Also, Bernoulli's principle can't give the reason for the 
yet another anomalous effect known as reverse Magnus 
effect 0. The right qualitative argument comes from the 
consideration of the boundary layer theory due to Lud- 
wig Prandtl[7| who hypothesized that for small viscos- 
ity, the viscous forces are negligible everywhere except 
close to the boundary where the no-slip condition has 
to be satisfied. The thickness of the boundary layer ap- 
proaches zero as the viscosity goes to zero. As the phe- 
nomenon of separation occurs on the line whose points 
are the singularities of the solutions of Prandtl's equa- 
tions, the boundary layer separates off the surface owing 



to the flow against the adverse pressure gradient and be- 
hind the sphere, there's the formation of the wake[H 
which is not axisymmetric for the rotating sphere; in this 
case there is upward displacement of the wake and ergo, 
a downward force on the sphere. 

III. THE MODEL EQUATION 

Anyway, it is no way possible to write a general 
equation for the swerving free-kicked football using the 
boundary layer theory. The saving feature is that the ac- 
tual force is very similar to the Magnus force calculated 
for the 1-D case. Consequently, after having explored 
the true reasons for the deflection, we adopt the Magnus 
force for the quantitative calculations. Hence, although 
a realistic free-kick is a 3-D affair, if the 1-D arguments 
as discussed in the beginning are applied to this case, 
we get the model equations. For that we assume for 
simplicity sans the loss of generality: 

1. The off-center kick causes the constant rotation of 
the football. 

2. The shape of the football doesn't get appreciably 
distorted. 

3. Initially, the ball is given no velocity component 
along y-axis and g acts along -z-axis. 

4. Ql.v at initial time t = and to is in x-z plane 
making an angle a with the z-axis. 

The first assumption needs a bit of justification that has 
been taken up in the appendix. By picking the right 
form for the drag force using empirical Cd vs Re curve 
for a sphere moving in a fluid at rest, one may write for 
10 3 < Re < 10 5 (which is typical of a free-kick) [io| : 

M'v = -k\v\v-Mg+ -ira 3 p(uj x v) (3) 
3 

where k is a constant and equation is subject to the initial 
condition at t = 0, (v x , v yi v z ) — (v X Q,v y o,v z o). In com- 
ponent form this can be broken down and rearranged as: 

v x = -A\f(v 2 x + v 2 + v\)v x - v y C cos a 

v y — — A\f(v x + Vy + v 2 z )v y + v x Ccosa + v z Cs'ma 

v' z = —Ay/(v%. + Vy + v 2 )v z — VyCsina — g (4) 

where A, C, a — tan _1 (u 2 o/w 2; o) are the constant param- 
eters. 

So, one has to use the set of equations as the model 
equation and fix A and C empirically which can even 
take care of reverse Magnus effect that can be effected 
by setting C as negative. The value of C obtained using 
Bernoulli's equation in within the order of magnitude of 
what is achieved using experiments [TlT| : so we shall stick 
to that for the sake of further discussion. We choose 
A = 0.027 and C = 0.16 (using Table 1). 
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Table 1: Values of constants (in SI units) 


Circumference of football (27ra) 


0.70 


TV yr £ £ j_1 11/71 T\ 

Mass ol loot ball (M) 


0.42 


Density 01 air [p) 


1.2 


Kinematic viscosity 01 air \y) 


14.6x10 


Typical velocity (v) 


on 
20 


Typical length scale (2a) 


U.zz 


Revnolrls nnmVier f /?p - — ( / ij)( c ?n)/i/\) 


~ 10 b 


Drag coefficient (Cd for 10 a < Re < 10 b ) 


-0.5 


Number of rotations on the ball (lo/2tt) 


2 


Acceleration due to gravity (g) 


9.8 


A of equation $4$ (A = 0.hpC D -na z /M) 


0.027 


C of equation Q (C = IMn^puj/M) 


0.16 




IV. THE FREEKICK IS NON-CHAOTIC 



v::: 





x-component of velocity (v x ) 



FIG. 2: Trajectories and fixed points: Going clockwise, 
Fig. 2a shows the locus of fixed points in 3-D, Fig. 2b (only 
v x — v z plane has been shown for clarity) is showing the col- 
lapse of various trajectories (dashed curves, starting form var- 
ious planes of the form v z = cv x marked in the figure) on the 
locus of fixed points (solid loop), and Fig. 2c is the blown up 
image of Fig. 2b near the locus of fixed points. 



Now, let's analyze the equations (|4]) purely as a set of 
non-linear dissipative (and hence, can neither have any 
quasi-periodic solution nor can have any totally repelling 
fixed point or orbit) coupled 3-dimensional differential 
equations. We set v = and solve the R.H.S to get 



FIG. 3: The negative Lyapunov exponent: The y- 

coordinate any arbitrary point on the curve is showing the 
average of the Lyapunov exponents calculated after every it- 
eration upto that point of the program [Ti|. Of course, the 
Lyapunov exponent is settling for a negative value and hence 
signifying that there should not be any chaos. 



following fixed points: 

v x = v* = —A^ 1 {gC 2 sinacosa} 

v y = v* = -A^ 1 {gHC sin a} 

v z = v* =+ A- 1 {g(H 2 + C 2 cos 2 i 



')} 



(5) 



where, A = H{H 2 + C 2 } and H is the real nega- 
tive solution of the equation H 6 + C 2 H 4 - A 2 g 2 H 2 - 
A 2 g 2 C 2 cos 2 a — which guarantees one value of H for 
each a. These are either spiral nodes or nodes for all a 
as can be obtained from the linear stability analysis (and 
extending the conclusions to this set of non-linear equa- 
tions, thanks to Hartman-Grobman theorem[10() which 
leads to following characteristic equation (A being the 
eigenvalues and £ = \/~(v* 



*2 , *2 

y 



+ 0): 
5A 2 £ 2 )\ + 2A 3 f 



A 3 + 4A£A 2 + (C 
AC 2 

+AC 2 £+ —(v* sin a- v* z cosa) 2 = (6) 

But note in the definition of (v* , v* , v* z ) there is a which 
being ta,n~ l {v Z Q/v x $) is a function of initial points and 
so, mathematically, for all a £ [0, 27r], we have a bounded 
locus of infinite number of fixed points (Fig. 2a) which 
being stable forces trajectories to converge on it at differ- 
ent points depending on which plane passing through the 
y-axis does the initial point of the corresponding trajec- 
tory lie (Fig. 2b). So, any initial point, however far, 
will evolve in time and collapse on the locus of fixed 
point, showing that the trajectories can't diverge, thanks 
to the bounded nature of the finite locus of fixed points. 
Hence, we may say that there is no chaos in the system 
of equations ((3]), for there's no exponential divergence of 
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the trajectories. This may be quantitatively justified by 
calculating, of course numerically, the largest Lyapunov 
exponent of the system of equations (j4|) and finding that 
the exponent is negative (Fig. 3), thereby confirming 
that there surely is no chaos in the velocity space of the 
football. One may ask, what happens if the ball slows 
down to below Re ~ 10 3 ? Can chaos come into picture? 
The answer probably is in negative since that is basically 
the Stake's law (with the possible Oseen's correction) 
dominated regime where the drag along each direction 
is linear in the velocity's component in that direction 
and of course, a set linear equations can't give rise to 
chaos. What remains unexamined is what happens pre- 
cisely when Re > 10 5 when the peculiar effect such as 
drag crisis Q and the compressibility effect are switched 
on. So what has been arrived at is in accordance with 
what we see in practice. A good free-kick specialist who 
has done enough homework can always almost send the 
ball where he wants to because the trajectory of the ball 
is grossly insensitive to the minor changes, natural even 
for the best of the players, in the initial velocities im- 
parted to the ball. 
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x-component of velocity (v x ) 



V. PECULIARITIES OF MODEL EQUATION 

Surprises continue to pour in for this simple system. 
Physically, a fixed point in the phase space means that 
if one starts from there, one remains there[l|. One more 
glimpse at the set of equations (0]) , and one can see that 
this long held intuitive idea for the fixed point has to 
be dispensed with at least for the dynamical systems for 
which there is a parameter depending on the initial con- 
ditions themselves. Any trajectory starting at a point 
lying on a plane v z — cv x (say) reaches a point on the 
locus of fixed points. Now, if one starts with the point on 
the locus of fixed points, the point is bound to be on some 
other plane v z = c'v x (say) and hence should evolve with 
time to reach another point on the locus (Fig. 4). So, as 
time tends to infinity a fixed point on the locus doesn't 
remain where it has been at t = 0, rather it moves to 
some other fixed point and hence in this case, it is the lo- 
cus of fixed points which maintains itself intact with time 
rather than the usually believed isolated fixed points. 



VI. CONCLUSION 

To conclude it may be highlighted that much to the de- 
light of mathematicians and others, it has been found in 
this brief paper that there exists, in the down-to-earth 
situation of ball-games, a class of dynamical systems 
(which, surprisingly, has been reported rather recently 
in the more-than-a-century-old literature on non-linear 
dynamics which pervades many fields viz. mathemat- 
ics, physics, chemistry, biology, economics etc.) where 
at least one of the parameters is dependent on initial 
points giving rise to the possibility that there may be 



FIG. 4: The peculiarity of the fixed points: The solid 
closed loop is the projection of the locus of fixed points on v x — 
v z plane. The curve on the right of point P is the trajectory 
for the evolution of initial point v x — 20, v z — 10 which has 
reached the loop at the point P that itself lies on the plane 
v z = ( — 18.8i>:,;)/0.681 (straight line passing through the point 
P) and hence, if is treated as an initial point, evolves as the 
broken curve on the left of point P to another point Q on the 
locus of fixed points. 



non-denumerably infinite number of fixed points. In the 
system discussed, the entire phase space got partitioned 
into a set of infinite number of basins of attraction. Also, 
the conventional intuitive definition of fixed point has 
been shattered because it has been seen here that an ini- 
tial condition on the fixed point can evolve as well on 
its own with time. Besides, it has been shown that a 
free-kicked football, which basically is such a dynami- 
cal system, cannot have chaos during its flight and hence 
the footballers may take free-kicks without bothering, too 
much precisely, about the initial velocity he should de- 
liver to the ball. This is in keeping with what happens 
in reality. 
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APPENDIX A: CONSTANCY OF THE 
ANGULAR VELOCITY 

Lets ponder upon the issue of constancy of the angular 
velocity. One might argue that the viscosity of the fluid 
might slow down the rotation of the football (or any other 
rotating sphere, to be general) as the time passes by and 
hence taking the angular velocity to be constant in the 
equation ^ is not justified. True, one knows that there 
is something called Kirchoff's law in fluids which gives 
the following expression[8[ for the instantaneous torque 
that slows a rotating sphere of radius a down: 



Sirrja Q 



(Al) 



where, rj is the dynamic viscosity and f2 is the instan- 
taneous angular velocity of the sphere. Obviously, if we 
denote the moment of inertia of the sphere by the symbol 
/, the torque-balance equation 



m = 



(A2) 



would yield the result that the angular velocity decays 
down exponentially fast. But, this conclusion hinges on 
the fact that the relation (|Aip requires the Reynolds 
number (Q,a 2 /v) to be much smaller than unity - a con- 
dition grossly violated in the case being considered in this 
paper. Using the table-1, we can see that 



n a 2 

v 



2(2tt) x 0.11 5 
14.6 x 10- 6 



10 4 



(A3) 



where we have taken initial value of Q i.e., to. Therefore, 
certainly relation (|A1|) is not valid for the case being con- 



sidered! 

So, what is the story here? For finding the analogous ex- 
pression for t in the case of high Reynolds number, lets 
resort to the powerful and yet simple method of dimen- 
sional analysis. The relevant parameters are p, a and Q. 
Simple dimensional analysis would then give 



5r>2 



-Npa Q 



(A4) 



Here, N is a numerical constant and the negative sign 
emphasises that this torque is trying to slow down the 
sphere's rotation. This expression might be given the 
physical interpretation of reaction on the sphere by the 
fluid due to the transfer of the angular momentum to 
the fluid by the sphere by dint of its rotation. Putting 
relation (|A4[) in the differential equation (|A2[) and solving 
with the initial condition: at t = 0, O = ui, one gets: 



n = 



i 



Npa s t 



(A5) 



Now, with some elementary knowledge in the non-linear 
dynamics, one knows that a fixed point is reached expo- 
nentially fast; relation (| A5|) suggests that the velocity of 
the rotational degree of freedom (i.e., the angular veloc- 
ity) slows down only as the inverse of time. Hence, one 
may conclude that as far as the behaviour of the trajec- 
tories, governed by the equations (UJ, to reach the fixed 
point is concerned, one can safely choose the angular ve- 
locity to be constant. 
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